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Abstract 

This paper proposes a self-optimized coverage coordination scheme for two-tier femtocell networks, 
in which a femtocell base station adjusts the transmit power based on the statistics of the signal and 
the interference power that is measured at a femtocell downlink. Furthermore, an analytic expression 
is derived for the coverage leakage probability that a femtocell coverage area leaks into an outdoor 
macrocell. The coverage analysis is verified by simulation, which shows that the proposed scheme 
provides sufficient indoor femtocell coverage and that the femtocell coverage does not leak into an 
outdoor macrocell. 

I. Introduction 

Femtocell technology has been emerging as a solution to the increase of both capacity 
and coverage while reducing both the capital expenditures and operating expenses of cellular 
networks. As femtocells share spectrum with macrocell networks, controlling the cross-tier 
interference between femto- and macrocells is need to be considered first in the enhancement of 
coverage and capacity. In addition, since a network operator may not be able to control femtocell 
locations, it is necessary for femtocells to sense the radio environment around them and carry 
out the self-configuration and self-optimization (HI- (Si of radio parameters from the moment 
they are set up by a consumer. 

Conventional dynamic cell sizing schemes, which adjust the transmit power of a base station 
(BS) [lUlISl or both the transmit power and antenna beam forming [|6l, have been developed 
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to improve the overall system capacity compared to that of a fixed cell sizing scheme. These 
approaches are not suitable for a macro/femto overlaid cell structure, where femtocell coverage 
must be controlled so it does not interfere with the outdoor macrocell. In order to achieve 
this goal, Claussen et al. JTl proposed a femtocell coverage coordination method that adjusts 
the femtocell pilot power, based on the number of handover events from outdoor passing, and 
the indoor users, which is robust against the varying size and shape of buildings. However, 
outdoor users may already experience inferior link quality during the time that a femtocell BS 
reduces its transmit power after recognizing the handover events of outdoor users. In particular, 
in a private access scenario that serves only registered users, the unauthorized users near a 
femtocell have a serious increase in the call drop rate or reduction of data rate. Moreover, the 
procedure that reconnects the rejected outdoor user to the macrocell may induce an additional 
handover, which causes a considerable amount of data transmission delay as well as packet loss 
in a packet switched cellular network with hard handovers, such as IEEE 802. 16e WiMAX and 
HSDPA imim. Therefore, this problem is severe for delay- and packet-loss-sensitive real-time 
applications such as Voice over Internet Protocol (VoIP). 

This paper proposes a coverage coordination scheme that is based on the statistics of the 
signal and interference power measured at a femtocell downlink (as opposed to the scheme 
based on handover events) to prevent any handover of outdoor users in advance. The proposed 
scheme comprises both the self-configuration and self-optimization of femtocell pilot transmit 
power. With a self-configuration function, a femtocell BS initiates its transmit power based on the 
measurement of interference from neighboring BSs in a manner that achieves a roughly constant 
cell coverage. The femtocell BS then performs a self-optimization function that continually 
adjusts the transmit power so that the femtocell coverage does not leak into an outdoor area 
while sufficiently covering an indoor femtocell area. 

II. Downlink Transmit Power Control 

This study considers a two-tier cellular network composed of overlaid macrocells and underlaid 
femtocells in which both cells use the same frequency channel. A femtocell BS is located at 
the center of a building with radius r^. Both BS and user are equipped with an omni-directional 
antenna. The femtocell BS creates cell coverage with radius r^, which is adaptively adjusted by 
the proposed transmit power control in order to correspond to the building area i.e., = r^,. The 
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transmit power control is composed from a two step procedure, where the femtocell BS initially 
self-configures its power and self-optimizes cell coverage by using transmit power control based 
on the measurements of radio environments. 

A. Initial Self-configuration 

A femtocell BS measures the average received power of pilot (over multiple frames to average 
out fast fading effect) from the neighboring macrocell and femtocell BSs on a neighboring BS 
list. The femtocell then chooses the strongest pilot power, Ib,max, among them. The femtocell 
BS configures its transmit pilot power such that the received pilot power from the femtocell 
BS and the strongest macrocell BS are identical on average at an initial cell radius of rj„j, i.e., 
Tf = Tint- The Appendix shows that Ib,max is nearly identical to the interference power from 
the strongest macrocell both measured by and averaged over the users located at the initial cell 
edge. Thus, the initial femtocell pilot power Pf^mi (dBm) is determined such that the femtocell 
BS power received at rj„j is equal to h^max, as follows: 

Pf,ini min (^Ib,max ~l~ L (Tj^j) , -Pmax) • (1) 

Here, Pmax and L are maximum femtocell pilot power and path loss, respectively. The initial 
self-configuration only provides the initial cell coverage of a femtocell, which is refined by the 
following self-optimized power control. 

B. Self-optimized Power Control 

A femtocell BS measures the level of other-cell interference /^(O) that is received from 
neighboring macrocells and femtocell BSs. The femtocell BS evaluates the received interference- 
plus-noise power Zu{0) = 10 logio(-^"(0) + where W is the thermal noise power. The 
femtocell BS collects the time-averaged received power Q^^^ in dBm, which is measured by 
each femtocell user during the ith iteration and is fed back to the currently linked femtocell 
BS. Based on the decision variable! F'^*) = Q^^^ — 2'„(0), where Q*'*^ is the averaged Q*-*-* over 
femtocell users, the transmit pilot power of the femtocell BS at the ith iteration is updated by 



min (pf + AP, Pmax) for r» < 



^(,+1) ^ I — y-f ^ , - maxy - ^ - th, 

' mm (pi''^ - AP, P^,,^ otherwise. 



-f '—^'^ ) ^ max 

'r''' gives a rough measure of spatially averaged carrier to interference-plus-noise ratio (CINR) over a femtocell area. 



Here, AP is the power control step in dB. pj*^^-* is determined by comparing F^*) with a threshold 
^th = To + Ta, where < Fa < TA.max- In order to make this power control scheme work 
properly, it is essential to set the threshold appropriately, that is, determining Tq and FA,max- 

1) Statistical Threshold Tq: Tq is obtained from the statistical characteristics of F^*^ Under 
the approximation that femtocell users are uniformly distributed over a circular femtocell with 
a radius of r/ at the zth iteration, the probability density function (PDF) of random variable 
D, which represents the distance between the femtocell BS and user, is given as /^((i) = 
2d/ (rj*^ —el), d G [eq, rj*^ ], where is the minimum D. Here, all radii and distances are in 
meters. Both the outdoor and indoor path loss in dB are modeled as L(D) = As + lOn log (D/dg), 
where n and Ag denote the path loss exponent and the path loss at a reference distance of 
ds = I m, respectively. For the outdoor-to-indoor path loss, the wall penetration loss Lp is added 
to L{D). When rj-*'' < Vb, the time-averaged received power of a femtocell user is given by 
Q{i) ^ pii) - L{D). Then, from the PDF of D, the PDF of g(^) is given as 



/Q„(g) = 10^^1nl0/(5n(r«'-5^)), qe [pf - L{rf), Pf - L^Sq) 

The expected value of Q^*^ is given by 

E [g(')] = 5n/ In 10 + Pf -A,- lQn{rf^ \ogrf - el log^o) / (rf ' - el). (4) 

Let denote Iu{P>) as the other-cell interference averaged over the users uniformly located on a 
circumference of radius of D centered at their femtocell BS. Since /„(0) ~ Iu{D) for D < O.Srffe j 
from the Appendix, /„(0) Iu{r^f); therefore, F^*^ is approximated as 

F« = g« - Z„(0) ^ E [Q«] - Z{rf), (5) 

where Z{r'^l^) = 10 log^g (j^u{r^f^) + . Furthermore, when ^ < r^, the average CINR of the 
user located at a femtocell edge is equal to a CINR threshold 7^/,, i.e., Q'-'H^/^) ~ ^('"/'*) = 1th- 
From this constraint, F*^*^ can be further approximated as follows: 

= 5n/ In 10 + 'yth - -Jy?^ log for rf < n, (6) 

where (a) follows from ([H) and the equation Q*-*-* = Pj^^ — L [D). This interestingly shows that 
F*^*) depends only on rj*\ and increases and converges to Fq = + 7t/i as rj-*'' increases (see 
the F(*) graph (i < 15) in Fig. [21). 



(3) 
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Fig. |2] describes an example of the coverage adaptation process that uses the power control 
scheme with Tfh = Tq. When rj*^ < {i < 15), F^*) is determined by ^ and it is less than Fq. 
Therefore, as the iteration index i increases, the transmit power Pj*-* increases, and the femtocell 
coverage rj"* extends to a building wall. The first time that r^^ is equal to rb {i = 15), F*^^^^ is 

(i) 

a little less then Fq, which leads to an increase in Pj- ' . However, contrary to the case of z < 15, 
the increase of transmit power no longer extends the femtocell coverage to the outdoors, i.e., 
rj*^^'' = r^j \ until the transmit power becomes large enough to overcome the additional path loss 
due to wall penetration so that femtocell coverage leaks into the outdoor region (see Case 2 in 
Fig. [3]). This constant cell coverage causes constant average interference-plus-noise power at the 
cell edge, i.e., Z{rf^^^) = Z{rf), thus F^^+^^-F^*) = E[Q(^+^)] -E[Q«] from ^. Additionally, 
/i+i) ^ ^« E[Q(^+i)] - ^[Q(^)] = - Pf from @, and P/+'^ - pf = ±AP 

from 0, i.e., F(*+i) - F^*) = ±AP. Therefore, F^^) is reformulated as 

I F(*) - AP otherwise 

According to this equation, F^^*^) = F^^^) + AP > Fq on the assumption that AP > Fq - F^^^) 
(more iterations will provide F*^*) > Fq when AP < Fq — F^^^^), and Pj*^ increases no more 
than Pf^^^- Thus, Tth, which is set to Fq, provides femtocell coverage that corresponds to the 
building area, i.e., rj = r^. It is important to note that this method is effective irrespective of 
the building's size, because Fq = + 'jth and F*^*) given by © or (|7]) are independent of r^. 

2) Maximum Additional Threshold FA,max-' From Fig. [21 we can observe that when the 
threshold Tth higher than Fq does not increase transmit power P|*^ up to the level at which 
femtocell coverage leaks into an outdoor area, a downlink CINR of femtocell is improved while 

(i) 

ry = Tb- From this observation we define FA,max as the maximum increase, which satisfies the 
condition rj'^ = r^, from the basic threshold Fq. 

rA,max is designed from the two casej^ of femtocell coverage described in Fig. [3l In Case 1, 
femtocell coverage is extended to a building wall, i.e., rj = rb, by using Tth = Tq. Increasing 
the transmit power of femtocell BS more than that of Case 1 results in Case 2, in which the 
femtocell coverage begin to leak into an outdoor region, i.e., rj = r^, + AD, where AD is a 

^In this section, several parameters (r^'\ Q'*', r^f \ and P^) defined in the previous section are classified into the parameters 
of Cases 1 and 2 using subscript numbers i and 2, and the iteration index ''^ is abbreviated for notational convenience 
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very small. Thus FA.max is given as 

TA.max = r2 — Tl (8) 

Note that while the femtocell BS increases the transmit power from Case 1 to Case 2, r/ remains 
equal to r^,, and Z(rf) is invariant. From ([5]), we then obtain 

Ta-Ti = E[Q2]-E[Qi] 

(a) 

= -P/,2 - Pf,l 

= Q2{n)-Qi{n) (9) 

where (a) follows from ^ and (b) follows from the equation Q{D) = Pj — L (D). In Case 1, 
as Fi = Tq and the CINR constraint Q{ri,) — Z{ri,) = '^th is preserved, the received pilot power 
of the femtocell at D = ri, is given by 

Qi{rb) = ith + z{n) (10) 

In Case 2, the boundary between the femtocell and macrocell is defined as the position where 
the received pilot powers of both cells are identical. Therefore, (52(^6 + ^-D) = /max(^b + '^D) 
at D = rf, + AD, where /max is the received pilot power from the strongest interfering BS. When 
this condition is satisfied, (52(^5) is given as 

Q2{rb) = Im^in) + 2Lp (11) 

Combining dH]), ©, (HO]), and ([lU), FA.max is given as 

TA.max = -^max(^6) + 2Lp - -fth " Z{rb) (12) 

In conclusion, Tth = Tq + Fa, where < Fa < FA,max, provides a higher downlink CINR 
than Tth = Fq due to the additional femtocell transmit power, while preserving the femtocell 
coverage that corresponds to the area of the building. 

III. Femtocell Coverage Analysis 

The statistical threshold Fq is derived from the expected value E [Q], but F, used for the 
power control, is estimated by using the sample mean Q at a femtocell BS. This results in a 
coverage leakage. Therefore, the coverage leakage probability that femtocell coverage leaks into 
the outdoor macrocell, Hk, is derived in this section. It is assumed that K femtocell users are 
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uniformly distributed in a building, and one of them is located at the boundary of the building, 
D = ri). The received pilot power averaged over K femtocell users is given by 

Q = (Ef=i'Qfc + Qk) /k= (Ef=i'Qfc + (Pf -As- 10nlog(r,))) / K, (13) 

where Qk and Qk is the received pilot power of the A;th femtocell user and the femtocell user 
located at D = rt, respectively. When r/ has approached rt, the T estimated in a femtocell BS 
is given asT = Q — Z(rb), and the femtocell BS increases its transmit power until T increases to 
Tth- If the additional femtocell transmit power, which is estimated to be Tth — T, is greater than 
TA.max, the femtocell coverage leaks into an outdoor area. Thus, Hk is defined and determined 
by using r^;, = To + Ta, Tq = E[Q] -Z{n) = E[Q]-Qk + -fth, and Z^{0) ^ Z{n), as follows: 

Hk = Pr [r^;, - r > rA,max] 
^ Pr [g < ^[g] + Ta - rA,max] 

= Pr[Q-gA^<ro-7t/. + rA-rA,max]. (14) 

Let the random variable X be defined as X = (Q — Qk) / K. Then, from the PDF of Q, the 
PDF of X is given as 

fx{x) = \2e-^'\ X G [0, 10nlog(rb/eo) / K] , (15) 

where Ai = iTlnlO/Sn, As = i^ln 10-10 '^""^Tn'""'''''' e-^Q^V(5n {rl - eg)) = r^iTln 10/(5n(r2- 
^o))- As A2 ~ Ai for rl ^ e^, the PDF of X is approximated to that of the exponential random 
variable with a parameter of Ai, i.e., fx{x) ~ Aic^'^^^. Next, let denote Ym as the sum of M 
independent, identically distributed random variables {X„i}m=i,--- ,m with a PDF identical to that 
of X: 

Ym = T.Z=iXm = T.Z=i {Qm - Qk) /K (16) 

The cumulative distribution function of Ym is then approximated to that of an Erlang random 
variable that was obtained by adding M independent exponential random variables with a 
parameter of Ai as ifTOl 

Fy,M ~ 1 - 5Z ^^r^' y ^ lOriMlog in/ So) /K] . (17) 

m=0 

From (O and ([M]), Yk-i = Qk - Qk, and ?/o = To - 1th + Pa - PA,max; (O is rewritten as 

Hk = Pi[YK-i<yo] 
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~ 1 - E^Zle-'^'^'i^iVor/ml , for > (18) 

and Hk = for yo < 0. 

In (fTSl ). e~^^y''{Xiyo)'^/m\ > because yo > 0. Thus, decreases as K increases, which 
demonstrates that a larger number of femtocell users improve the performance of the proposed 
scheme. Additionally, from the Taylor series for the exponential function = X]m=o^'"/^'' 
the asymptotic behavior of Hk at a very large K is given as 

lim Hk = 1- e-^^^« y ^^i^ = 0. (19) 

This indicates that the proposed scheme is asymptotically optimal in terms of coverage leakage 
probability. 

IV. Performance Evaluation 

A. Simulation assumptions and performance metrics 

With the system parameters given in Table I, a Monte Carlo simulation approach is used to 
evaluate the coverage leakage probability Hk, the average leakage distance and the average 
indoor coverage ^. A macrocell has a layout of 19 hexagonal cells arranged in a hexagonal 
lattice with two rings of cells surrounding the center cell. A target femtocell is located at a point 
with distance di, = 400 m from the macrocell BS in the center cell, and 50 interfering femtocells 
with a fixed pilot power given from ([T]) are uniformly distributed within the center macrocell with 
a radius of = 580 m. For each simulation repetition in a Monte Carlo simulation with 5,000 
trials, the indoor and outdoor users are uniformly distributed in the building and the outdoor 
area from 10 m of the building wall, respectively. After the geometrical configuration of the BS 
and users, the femtocell pilot power is initiated from ([T]) and is optimized by using ^ until it 
converges. The femtocell coverage is then evaluated so that the received pilot power of femtocell 
is larger than that of the macrocell. Hk is obtained by dividing the number of events where 
femtocell coverage leaks outdoors by the total number of trials, n is evaluated by averaging the 
distance between the leaked femtocell edge and building wall over the total number of trials. ^ 
is calculated as the ratio of average femtocell coverage to the building's area. Hk and measure 
femtocell coverage that leaks into the outdoor area. Additionally, ^ measures the performance 
of providing sufficient indoor femtocell coverage. 
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B. Simulation results 

Fig. |4] shows Hk versus Fa for K=5, 10, 20, 40, and infinity wlien Lp = 10 and 3 dB. The 
analytic curves are very close to the simulated curves. Hk increases as Fa increases due to 
the fact that the higher Tth that is induced by an increase in Fa causes a rise in the transmit 
power of a femtocell BS. A larger wall penetration loss Lp reduces the Hk due to the increasing 
value of Fa, max as shown in (fT2l) . Moreover, Fa guaranteeing some level of Hk increases as 
Lp becomes larger, which results in the higher CINR of a femtocell downlink. Therefore, the 
proposed scheme is more effective for a building with a higher wall penetration loss. Fig. |4] 
also shows the impact of K on Hk- A larger K increases the probability that the sample mean 
Q becomes close to the expected value of E[Q], i.e., the approximation error of Q decreases, 
which results in less Hk where Hk is zero for infinite K. A description of this result is also 
in the last paragraph of Section |lVl At least 40 users (less users are probably deployed in a 
femtocell except for the enterprise scenario), should be deployed in order to achieve an Hk value 
of 5 % when Lp = 3 dB. However, according to available literature [fTT ] |[[T2 l| . the probability 
that Lp = 3 is low in actual applications. Thus, the proposed scheme remains preferable. The 
analysis and simulation of Hk as shown in Fig. |4] are performed by considering no shadowing 
and the perfect estimation of path loss exponent n. Although this assumption is unrealistic, it 
provides good insight regarding the proposed algorithm's performance according to many factors 
of Tth, Lp, n, and K. 

In a real femtocell scenario we consider shadowing, path loss exponent estimation error, a small 
number of users, and rectangular- shaped building with several non-symmetric wall penetration 
loss. Although shadowing is included in the path loss model, Fq does not change because the 
shadowing averaged over an indoor area is zero, i.e., E[Q] does not vary. On the other hand, more 
uneven cell coverage, due to both shadowing and several non-symmetric wall penetration loss, 
highly increases Hk, but despite of this higher Hk, the leakage area can be small so that few 
outdoor users are linked to a femtocell BS. Thus, average leakage distance Q and average indoor 
coverage \1' are investigated for the real scenario with parameters in Table I, and the results are 
shown in Fig. |5l The path loss exponent error is considered by using and n^, which denotes 
the path loss exponent n used for determining the Fq and calculating the received power in a 
real link, respectively. A higher Fa leads to a rise in the transmit power of the femtocell BS, 
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which increases and Thus, Fa is adaptively determined according to both the maximum 
permissible n and minimum achievable ^. From Fq = + 'jth, > increase Tth, which 
gives the same impact on the performance as Fa increases under the condition = Ur- On the 
contrary, < Ur leads to opposite results. Thus, a higher, or a lower, n,, is recommended to 
increase or decrease Vl, respectively. The simulation results indicate that the proposed algorithm 
archives less than 5 m and ^ more than 0.9, which is a feasible level of performance for 
the realistic scenario. Additionally, this scheme requires additional uplink overhead for reporting 
the average received pilot power. As the power is averaged over multiple frames to remove fast 
fading, long-term reporting sufficiently supports the amount of feedback information, i.e., the 
overhead is not so considerable as to make implementation impossible. 



This paper proposes a novel coverage coordination scheme based on a self-configuration and 
self-optimization of transmit power. An analytic expression for the coverage leakage probability 
of the femtocell is derived and verified by simulations. The simulation results show that, by using 
the proposed scheme, femtocells provide sufficient indoor coverage and low coverage leakage 
to outdoor area. In conclusion, the proposed scheme can make femtocell coverage correspond 
to the building's area without knowing about the area of the building. Further research needs to 
improve the robustness against the small number of users as well as to investigate the effect of 
mobility of users. 



Fig. \T\ shows the geometric configuration for calculating receiving power of the other-cell 
interference at a femtocell BS and at its users. It is assumed that all M interfering BSs use 
the same transmit power level in dBm and that the wall penetration loss Lp is constant 
irrespective of 9. The other-cell interference measured by and averaged over the users, which 
are assumed to be located uniformly on a circumference with a radius of D, is given by 



V. Conclusion 



Appendix 




(20) 
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where the subscript hn represents a linear value, and 2-^1 [■] is Gauss's hypergeometric function 
[fT3l . Iu{D) can be obtain in a further simple form when n = 2k, as follows: 

2fc — 1 

UD) = J2{D + d,,r'''-\D^ + dir'/ (i - ' , (21) 

AsM-'^pM^i ^ [J-^ + db,i) J 

where is a integer [fT3ll . On the other hand, the other-cell interference at the femtocell BS is 
given by 

= ^^Xt^'^-^l (22) 
Let \J and B denote the summation part of (|2TI) and (|22|) . respectively. If n = 2, 4 and (ift j = 
2/}, 3-D, and AD for all z, f/ and i? are given as follows: 



\J 



0.333M 


,'0.185M 


1)2 , 




0.125M 


/0.019M 






0.067M 


/' 0.005Af 


£)2 5 





) 

) > ^ 
) 



0.25M 


,'0.063Af 






O.lllAf 


( 0.012A/ 






0.063 A/ 


/' 0.004Af 


£)2 , 


I 1)4 



= 2D 

db,i = 3D forn = 2,(or4) (23) 

= 



The similar values of U and 5 shows that Iu{D) ^ /n(0) for the realistic conditions of a path 
loss exponent n larger than 2 and cife^j > 2D. 
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TABLE I 

System parameters for the simulation and analysis 



Parameter 


Symbol 


Value 


Macrocell radius 




580 m 


Macro-to-femto BS distance 


db 


400 m 


Building sliape 




Theory verification (Fig. 4): Circle 
with a radius of rb = 20m, 
Real scenario (Fig. 5): 
20m X 15m rectangle 


Wall penetration loss 
(Percent value represents 
wall-length ratio) 




Theory verification (Fig. 4):Lp = 3, lOdB 
Real scenario (Fig. 5): 
15dB(35%), 10dB(30%), 7dB(20%), 2dB(15%) 


Initial femtocell radius 




15 m 


Minimum distance 
between a femtocell BS and a user 


£0 


1 m 


Path loss at 1 m 


As 


37 dB 


Power control step 


AP 


0.25 dB 


Maximimi transmit power of a BS 




23 dBm (femto), 43 dBm (macro) 


Thermal noise power 


w 


-96.8 dBm 


Path loss exponent 


n 


3 (femtocell), 4 (macrocell) 


CINR threshold 




-2.6 dB 


Statistical threshold 


To 


3.91 


Maximum additional threshold 


rA,max 


1.52 dB (Lp=3), 9.8 dB (Lp=10) 
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Building wall 



Fig. 1. Geometric configuration for calculating the othercell interference at femtocell BS and users 
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Fig. 2. A change of transmit power, femtocell coverage, and F'*', using the proposed power control with Vtu = To and AP=2 
dB. 



15 




Fig. 3. Wall penetration loss Lp enables a femtocell BS to utilize the additional threshold Fa, < Fa < Fa, max (Case 1: 
r/ = rb. Case 2: rf — n + AD). 




Fig. 4. Hk versus Fa, comparing the simulation and analysis results for Lp = 10 and 3 dB. 
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Fig. 5. f2 and ^ versus Fa with shadowing , path loss exponent error (rie = 3), small number of users {K = 2), and several 
non-symmetric wall penetration loss. 



